Indications of Spin-Charge Separation at Short Distance and Stripe Formation 
in the Extended t-J Model on Ladders and Planes 
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The recently discussed tendency of holes to generate nontrivial spin environments in the extended 
two-dimensional t-J model (G. Martins, R. Eder, and E. Dagotto, Phys. Rev. B60, R3716 (1999)) 
is here investigated using computational techniques applied to ladders with several number of legs. 
This tendency is studied also with the help of analytic spin-polaron approaches directly in two 
dimensions. Our main result is that the presence of robust antiferromagnetic correlations between 
spins located at both sides of a hole either along the x or y axis, observed before numerically on square 
clusters, is also found using ladders, as well as applying techniques based on a string-basis expansion. 
This so-called "across-the-hole" nontrivial structure exists even in the two-leg spin-gapped ladder 
system, and leads to an effective reduction in dimensionality and spin-charge separation at short- 
distances, with a concomitant drastic reduction in the quasiparticle (QP) weight Z. In general, it 
appears that holes tend to induce one-dimensional-like spin arrangements to improve their mobility. 
Using ladders it is also shown that the very small J/t~0.1 regime of the standard t-J model may be 
more realistic than anticipated in previous investigations, since such regime shares several properties 
with those found in the extended model at realistic couplings. Another goal of the present article is 
to provide additional information on the recently discussed tendencies to stripe formation and spin 
incommensurability reported for the extended t-J model. These tendencies are here illustrated with 
several examples. 

PACS numbers: 74.20.-z, 74.20.Mn, 75.25.Dw 
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I. INTRODUCTION 

The study of high temperature superconductors con- 
tinues attracting the attention of the Condensed Mat- 
ter community. In recent years, much effort has been 
devoted to the understanding of the spin incommensu- 
rabihty that appears in neutron scattering experiments 
for some of these compounds |l| . Tranquada's interpre- 
tation of the experiments is based on "stripes" where 
charge is confined to one-dimensional (ID) paths in the 
crystal, with an average hole charge n;i=0.5, namely one 
hole for every two sites along the stripe [0. The stripe 
interpretation appears robust in the one-layer material 
La2-xSrxCu04, although it is still controversial in the 
bilayer YBa2Cu306+5 [0. Experimental results compat- 
ible with metallic stripes have also been reported using 
other techniques H. On the theory front, early stud- 
ies discussed the presence of stripes in Hubbard and t-J 
models, although typically with density n;i=1.0 [0jq|. Re- 
cently, the existence of n^, < 1 stripes was discussed using 
computational techniques directly in the t-J model at in- 
termediate coupling J/t, without the need of long-range 
Coulomb interactions M. Sometimes these stripes are 
described as a condensation of d-wave pairs |Q. Stripes 
appear also in Monte Carlo studies of the spin-fermion 
model P]. However, the origin of such complex spin- 
charge arrangement is still under much discussion, and 



even its stabilization remains controversial in the stan- 
dard t-J model ip. 

In parallel with the developments in neutron scat- 
tering experiments, in recent years photoemission tech- 
niques have provided a plethora of information about the 
cuprates. In particular, the experimental study of the 
one-hole spectral function of the parent insulator com- 
pound gave us useful information to judge the quality 
of existing models for these materials [Q. Based on 
the experimental one-hole dispersion it has been convinc- 
ingly shown that the t-J model is not enough to address 
the cuprates but extra terms must be added in order 
to reproduce the photoemission results for the insulator 
pd| . These terms appear in the form of extra hole hop- 
pings, regulated in intensity by hopping amplitudes usu- 
ally denoted by t' and t". The bare value of these extra 
amplitudes is small compared with the nearest-neighbor 
hopping t, but its influence is substantial since they link 
sites belonging to the same sublattice and they are not so 
heavily renormalized to smaller values as it occurs with 
t. The origin of these extra hoppings was discussed be- 
fore ||ll|] and they appear naturally in mappings from the 
three-band model for cuprates to a one-band Hamilto- 
nian, and also in band-structure calculations. In fact, 
what is unnatural is to assume that the standard t-J 
model, with t'=t"=0.0, is valid for the cuprates with- 
out corrections. Unlike in gauge theories, there are no 



renormalizability or symmetry arguments at work pre- 
venting the introduction of nonzero extra hoppings in 
models for copper oxides. For these reasons it is impor- 
tant to study the extended t-J model in detail and, in 
particular, whether stripes occur in its ground state. 

Recent studies of the extended t-J model have provided 
useful and sometimes surprising information |L2| , pj3[ . For 
instance, it was observed that the one-hole QP weight 
Z is considerably reduced by the addition of t' and t" 
at fixed J/t. This result is particularly dramatic at mo- 
menta (0,7r) and (tt, 0), and its origin lies in the gener- 
ation of across-the-hole antiferromagnetic (AF) correla- 
tions, namely in the reference frame of the mobile hole 
the two spins immediately next to it along one axis are 
aligned antiferromagnetically [[l2| . This tendency is the 
opposite as expected from a vacancy in an antiferromag- 
net where those two spins should be ferromagnetically 
aligned, since they are in the same sublatticc. Such a cu- 
rious result was also observed in previous DMRG studies 
of ladders and planes using the standard t-J model by 
White and Scalapino jlj] . Its importance and origin was 
recently discussed by Martins et al. [|l2| where its pres- 
ence was conjectured to be caused by tendencies in the 
system to spin-charge separation at short distances. In 
other words, the kinetic energy of the hole with the ex- 
tra mobility induced by t' and t" dominates most of the 
physics of the one-hole problem and it tries to generate 
in its vicinity an environment that allows for the hole 
to move without fighting against the spin background. 
In one-dimension this is easy to set up, and indeed an 
antiferromagnetic across-the-hole environment is stabi- 
lized on chains leading to spin-charge separation [[l5[ . 
In two-dimensions the hole tries to create such an en- 
vironment in its vicinity but it cannot produce a total 
spin-charge separation since frustration is induced when 
antiferromagnetic across-the-hole interactions are gener- 
ated along both axes. It is expected that for realistic 
values of the couplings, the spin-charge separation will 
only be local (i.e. at short distances ) [[l2[ . This result, 
obtained using small square clusters |12|], is clear in the 
presence of nonzero t' and t" of the proper sign and mag- 
nitude, but it also occurs in the standard t-J model in 
the regime of very small J/t, which remains mostly un- 
explored. This interesting observation redefines the rele- 
vant value of the parameter J/t from the number usually 
used in most studies (0.3-0.4) to smaller values close to 
0.1. 

Also motivated by the results in Rcf. ||l^ , recently the 
analysis of many holes at couplings that are expected to 
lead to robust across-the-hole correlations was reported 
by Martins et al. [|l3[. Indications of stripes were ob- 
tained in this case, similarly as those previously found in 
Ref. J6J, rationalized as the natural tendency of locally 
spin-charge separated states to avoid the spin frustra- 
tion caused by the nontrivial spin environment that each 
hole generates. The hole density in these stripes is in 
good agreement with experiments, and qualitatively the 
results resemble those in the "holons in a row" picture 



of Zaanen Q . Tight hole bound-states are not needed to 
generate stripes. 

These observations suggest that either the extended t- 
J model or the standard t-J model at very small coupling 
must be analyzed in detail if the goal is to reproduce the 
physics of the striped regime of the cuprates. In addition, 
it is important to confirm the generality of previous re- 
sults |12| by studying other systems with similar physics 
but where computational studies can be carried out on 
larger clusters. For this purpose, here ladders are inves- 
tigated using the extended and standard t-J models in 
the coupling regime where in Rcf. ||l^ the robust across- 
the-hole feature was observed on small square clusters. 
Our main result is that these robust correlations are also 
very clearly identified on ladders, even with only two legs, 
and, thus, its presence is more general than naively antic- 
ipated. Based on our results it is clear that holes on lad- 
ders tend to form onc-dimcnsional (ID) like environments 
in their vicinity to improve their kinetic energy leading to 
an effective reduction in dimensionality. Spin incommen- 
surability is generated by this procedure, a novel result 
in two-leg ladder systems to the best of our knowledge, 
although its presence in two dimensions has been known 
for some time [EgI. 

The present paper has a second goal which is the de- 
tailed analysis of stripe formation upon doping of ladders 
with many holes. For this purpose here the recent results 
of Rcf. ||l^] where stripes were observed are extended to 
couplings and parameters not reported before, to illus- 
trate and confirm the tendencies toward stripe formation 
in this context. The proposed rationalization for this 
stripe stabilization was already discussed |1^ and it is 
based on the sharing of the locally spin-charge separated 
spin environment created by each individual hole, thus 
avoiding spin frustrating effects. 

The model used here is the extended t-J model defined 



as 
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where tjm is t for nearest-neighbors (NN), t' for next 
NN, and t" for next next NN sites, and zero otherwise. 
The scale of energy will be t=l unless otherwise spec- 
ified. The rest of the notation is standard. Compari- 
son with PES experiments [Q showed that t'=-0.35 and 
t"=0.25 are relevant to explain PES data [|l|,|2|. To 
simplify our studies the ratio t'/t" will be fixed to -1.4 in 
most of our analysis, although other ratios will be used in 
some cases. As numerical techniques, the Density Matrix 
Renormalization Group (DMRG) ||q,p7[|. Exact Diagonal- 
ization |16|| , and an algorithm using a small fraction of 
the ladder rung-basis (optimized reduced-basis approx- 
imation, or ORBA 1 18 1) are here used. The paper is 



organized by the number of legs of the ladders consid- 
ered, from two to six. Analytical results are discussed 
after the numerical methods, illustrating the appearance 
of tendencies toward across-the-hole antiferromagnetism 



in two dimensional systems. In the last section it is con- 
cluded that the mostly unexplored extended t-J model 
contains interesting physics related to the cuprates, that 
deserves further work. 



II. TWO-LEG LADDERS: 
A. One Hole: 

The present computational analysis of ladder systems 
starts with the two-leg ladder case. Previous literature 
has shown that in the absence of t' and t" hoppings, and 
at intermediate values of J, such as 0.4, one hole behaves 
as expected from a carrier in an antiferromagnet even if 
the ladder spin background has only short-range mag- 
netic order p9[ |. In other words, a hole creates a spin 
distortion in its vicinity (spin polaron) and the one-hole 
ground state has a finite QP weight Z. Two of these spin 
polarons bind into a hole pair that leads to supercon- 
ducting correlations upon further doping |19|| . However, 
in this paper our goal is to analyze a region of parameter 
space not studied before, to the best of our knowledge, 
where the hole is expected to be substantially more mo- 
bile. For this purpose J is reduced to values between 0.1 
and 0.2, still within the domain considered realistic in 
studies of the cuprates, and t', t" are made nonzero and 
of the sign and magnitude as suggested by photocmis- 
sion experiments. At such "small" values of J/t, the spin 
background no longer is expected to fully control the be- 
havior of the hole but the spins have to arrange in such 
a way that the hole kinetic energy is optimized (but still 
without leading to a ferromagnetic state that may be the 
best configuration at J/t^O.O). 

Consider first Exact Diagonalization results in the one- 
hole sector. In Fig. la, the spin-spin correlations are 
shown for the case where the hole is projected on the 
(arbitrary) site indicated, from the full one-hole lowest- 
energy state with momentum k=(7r, 0), working at J=0.2, 
t'=-0.35 and t"=0.25. It can be observed that at a dis- 
tance of about three lattice spacings from the hole the 
spin correlations are similar to the ones expected for an 
undoped ladder, with robust antiferromagnetism along 
both the rungs and legs. However, in the vicinity of the 
hole the spin correlations are substantially altered. In 
this case the spins belonging to the same rungs are no 
longer strongly antiferromagnetically correlated and the 
system in the vicinity of the hole resembles a pair of 
weakly-coupled portions of a chain, one per leg. In this 
respect it appears that an "effective" transition from a 
two-leg ladder to ID chains has occurred near the hole, 
a surprising result. Holes seem to optimize their energy 
by creating ID environments in its vicinity. Note the 
presence of robust across-the-hole spin correlation on the 
upper leg of Fig. la, which also appears in the exact solu- 
tion of the ID Hubbard model due to spin-charge separa- 
tion. Fig. la suggests that the spin associated to the hole 



is spread uniformly in the ID-like segment dynamically 
generated, and in this respect spin-charge separation oc- 
curs locally, as anticipated in the Introduction and as 
it was reported in Ref. [Q using small square clusters. 
Note in Fig. la that the spin belonging to the same rung 
as the hole appears as "free" and there is a strong AF 
bond at distance of two lattice spacings across it, as it 
occurs across-the-hole. 
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FIG. 1. Spin-spin correlations in the state with the lowest 
energy of one hole with momentum k, at the couplings de- 
scribed below. The data were obtained on the two-leg ladder 
using the Lanczos method on 2x8 (a-d) and 2x 12 (e) clusters, 
with periodic boundary conditions (PBC) along the leg direc- 
tion. Shown are results when the hole is projected from the 
state under consideration at the site indicated by an open cir- 
cle. The dark lines between sites represent the absolute value 
of the (antiferromagnetic) (Si ■ Sj) spin correlation between 
the spins located at those sites (scale shown in (a)). Dashed 
lines indicate weak ferromagnetic correlations. Results at dis- 
tance of two lattice spacings are shown only near the hole. 
Note the presence of strong across-the-hole antiferromagnetic 
bonds in all the cases, (a) corresponds to J=0.2, t'=-0.35, 
t"=0.25 and k=(7r,0). Far from the hole the spin correlations 
are as in undoped two-leg ladders, while near the hole they 
are much distorted, and they are robust mainly along the legs, 
(b) is the same as (a) but for t'=t"=0.0. (c) Results at J=0.1, 
t'=-0.35, t"=0.25 and k=(7r,0). Note the appearance of some 
weak ferromagnetic rungs at this small value of J. (d) Same 
as (a) but for k=(7r/2,0). (e) corresponds to J=0.2, t'=-0.35, 
t"=0.25 andk=(7r,7r). 

The comparison of Figs. la and lb, the latter obtained 
without extra hoppings, shows that t' and t" are impor- 
tant at J=0.2 to regulate the size of the ID environment 
around the hole: as t', t" grow in magnitude, the size of 
the ID-like region also grows. This same effect is pro- 
duced reducing J at fixed t',t" as shown in Fig.lc: here 



the distortion is larger than at J=0.2 and even weak ferro- 
magnetic Hnks are generated between the chains, an un- 
expected result quite different from the well-established 
physics of undoped ladders at intermediate J. The re- 
sults here are also weakly dependent on the momentum: 
in Fig. Id, the case of k=(7r/2, 0) (momentum of the over- 
all ground-state of one hole) presents a distortion similar 
in size as at k=(7r,0), although the AF long-bond in the 
leg opposite to where the hole is projected (Fig. la) is 
no longer present. However, the long-bond across- the- 
hole is still robust. Fig.le is another case illustrating the 
momentum dependence in the problem, corresponding 
this time to k=(7r, tt). Overall, it can be concluded that 
the lowest one-hole energy states for all momenta stud- 
ied here, and in large regions of parameter space, present 
very similar spin arrangements. The across-the-hole fea- 
ture emphasized in Ref. |12|] is very robust in all cases 
and in this sense there is a clear 7r-shift across- the-hole 
even in the simple case of a two-leg ladder, similarly as it 
occurs in 2D systems. The dynamically induced ID-like 
regions near the carrier are also clear in our studies. 
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FIG. 2. One-hole spectral function A(k,(j;) obtained ex- 
actly from the undoped ground state of a 2x12 cluster work- 
ing at J=0.2, and t'=-1.4t" fixed, parametric with t'. The 
arrows indicate the position of the first pole in the spectra 
and the numbers next to them are the weights Z (normal- 
ized such that 0<Z<1). (a) corresponds to k=(7r,0), (b) to 
k=(7r/2,0), and (c) to k=(7r,7r). 

The generation of a quasi- ID structure in the vicinity 
of the hole and the expected spread of the hole's spin- 
1/2 over several lattice spacings (spin-charge separation 
at short distances) causes a drastic reduction in the QP 
weight Z. In Fig. 2, the one-hole spectral function A(k,a;) 
is presented at three different momenta. The results are 
shown as it is usual in experimental photoemission liter- 
ature with the lower energy states appearing near cl;=0. 



reference energy which is located at an arbitrary position 
in this half-filled case, and the rest of the states running 
to the left. Data for J=0.2 and several values of t', at 
fixed t'/t"=-1.4, are presented, together with the weight 
Z and position of the first state (arrows). Fig. 2a shows 
that Z at this value of J and k=(7r, 0) is already small 
even in the absence of extra hoppings, result similar to 
those observed in previous investigations of the t- J model 
on square lattices |16[. However, it is remarkable the 
rapid reduction of Z with increasing |t'|, to values that 
in practice are basically negligible by the time these ex- 
tra hoppings reach its realistic values. This is correlated 
with the appearance of the complex quasi-lD structure 
discussed in Fig.l [E0| A similar phenomenon but even 
more pronounced exists for k=(7r, tt) (Fig. 2c), where Z 
is negligible even without extra hoppings. On the other 
hand, the results at k=(7r/2,0) (Fig. 2b), while still cor- 
responding to small values of Z of the order of 10% of the 
maximum, do not show such a dramatic reduction with 
increasing t'. In this case the spin appears still located in 
the vicinity of the hole, although it is spread over several 
sites. Z becomes negligible at this momentum only at 
values of |t'| larger than shown in the figure. 



0.20 



0.15 



(a). Z^- 



■ ■ ■ 



° 1.5 
« 1.0 
N 0.5 






~V-V' 



N 



0.10 



0.05 



0.00 



0.0 0.2 0.4 0.6 0.8 1.0 

,n — D |f/t| 
(7i/2,0)\^^ 

d" 



6^371/4,0) 

^°\ 

■ I ¥ 




0.0 0.2 0.4 0.6 0.8 1.00.0 0.2 0.4 0.6 0.8 1.0 

It'/ti itvti 

FIG. 3. (a) Weight Z of the lowest energy state in the 
A(k,tj) one- hole spectral function of exactly solved 2x12 clus- 
ters at J=0.2, t'=-1.4t", as a function of |t'/t|. Z is normal- 
ized such that 0<Z<1. Shown are results at the momenta 
indicated. In the inset the same information is given but for 
the nearly negligible lowest-pole weight of k=(7r, tt). (b) Lines 
of constant Z weight (value indicated) in the plane J/t-|t'/t|, 
with t'=-1.4t" fixed, corresponding to the lowest-energy pole 
in the k=(7r,0) subspace. Note the drastic reduction in Z as 
J is reduced from 0.4 and/or |t'| is increased from 0.0. The 
results are interpolations using a fine grid of points obtained 
from the 2x8 cluster solved exactly. 

The results for the QP weight Z are more explicitly 
illustrated in Fig. 3a where its |t'| dependence is shown 
for several momenta at fixed J=0.2. Clearly for all mo- 
menta the presence of t' eventually leads to a drastic re- 
duction of the Z value, although for k=(7r, 0) it happens 



rapidly with increasing t' while at k=(7r/2,0) it needs a 
t' larger than the nearest-neighbors hopping t. In Fig. 3b 
the lines of constant weight Z corresponding to k=(7r, 0) 
are shown. It is clear that Z at this momentum rapidly 
reduces its value, both at fixed J increasing |t'| or for fixed 
hoppings decreasing J. Although a careful finite size scal- 
ing is difficult, previous experience with two-dimensional 
clusters suggest that Z tends to decrease as the lattice 
size grows. As a consequence it is expected that the val- 
ues reported in Figs. 2 and 3 will actually be even smaller 
for a very long two-leg ladder. Note also that the lines 
of constant Z suggest that the physics of, say, interme- 
diate J and finite t' may be smoothly connected to that 
of small J and zero t', quite similarly as it happens on 
square clusters |12|| . Then, it may occur that analyzing 
the regime of abnormally small J (~0.1) of the t-J model 
may effectively account for the presence of extra hop- 
pings, as was conjectured in Ref. |12|. This assumption 
is interesting due to the anticipated unusual properties 
of the small J/t region of the t-J model, and it avoids the 
somewhat aesthetically unpleasant use of extra hoppings 
in the present theoretical studies. 



B. Several holes: 

The previous subsection showed that the two-leg lad- 
der system presents exotic behavior in the one-hole sector 
when the regime of small J is investigated or alternatively 
when at a fixed J the hoppings t'-t" are switched on. It 
is important to explore what occurs for more holes in the 
problem, and Fig. 4 contains exact results for the case of 
two holes on an 2x8 ladder. Data is shown for the case 
when the two holes are projected to its most probable 
location in the ground-state, which in the cases stud- 
ied coincide with the maximum distance among them 
allowed in the cluster. The spin structures in Fig. 4 near 
the holes clearly resemble those found around individ- 
ual holes in the previous subsection. In particular the 
across-the-hole AF correlation is very prominent in all 
cases shown in Fig. 4. This structure naturally leads to 
spin incommensurability, as observed in the spin struc- 
ture factor presented in Fig. 5a ||2^]. It is interesting to 
observe that the coupling between chains can be weakly 
ferromagnetic in some cases (Fig. 4a) , weakly AF in oth- 
ers (Fig. 4b), or even virtually negligible (Fig. 4c) depend- 
ing on the actual couplings used. Overall it is clear that 
the legs are approximately decorrelated upon hole doping 
at the couplings analyzed here, and each leg behaves as a 
ID-chain described by the same original t-J model. Each 
individual leg resembles the expected result for the ID 
Hubbard model at large U [^, as discussed before for 
2D clusters [12]. It is remarkable that very recently nu- 
clear spin relaxation results for two-leg ladder materials 
have also reported the decoupling of the two legs at high 
temperature |2^. The analysis of the relation between 
our results and those of Ref. 01 deserves further work. 
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FIG. 4. Spin-spin correlations in the state with the low- 
est energy of two holes, at the couplings described below. 
The data were obtained exactly on the two-leg ladder us- 
ing the Lanczos method on 2x8 clusters with PBC along the 
legs. Shown are results when the holes are projected from 
the state under consideration at the sites indicated by open 
circles. These hole positions are the most probable in the 
two-holes ground-state. The dark lines between sites rep- 
resent the strength of the antiferromagnetic spin correlation 
between the spins located at those sites (scale shown in (a)). 
Dashed lines indicate weak ferromagnetic correlations, and 
the missing vertical correlations indicate links with virtually 
negligible antiferromagnetic spin correlations. Results at dis- 
tance of two lattice spacings are shown only near the hole. 
Note the presence of strong across-the-hole antiferromagnetic 
bonds in all the cases, as in the one-hole problem. Note also 
the clear tendency of spins in opposite legs to be nearly decor- 
related, (a) corresponds to J=0.2, t'=-0.35, and t"=0.25. 
(b) corresponds to J=0.2, and t'=t"=0.0. (c) corresponds to 
J=0.4, t'=-0.35, and t"=0.25. 

In the two-hole ground-state the expectation value of 
the number operator for a particular momentum is shown 
in Fig. 5b. Clearly the momenta k^ dominating in this 
state is not only 7r/2, the momentum of the lowest energy 
state of the one-hole sector, but actually 37r/4 and tt are 
also important as well as the ky=7r sector. Then, there 
is no indication of "hole-pockets" in the doped system, 
and momenta over a wide range contribute to the two- 
hole ground-state. This suggests that the one-hole states 
with the largest weight in the two-hole ground-state have 
the across-the-hole structure and local spin-charge sepa- 
ration. 
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FIG. 5. (a) Spin structure factor S(kx,ky=7r) vs kx for the 
couplings J and t' indicated (with t'=-1.4t" fixed), using the 
two holes ground-state of the 2x8 cluster. Clear tenden- 
cies toward incommensurate correlations are observed in our 
study, as exemplified here, either at small J without t'-t" 
hoppings, or increasing the values of the latter at fixed J. (b) 
Exact n(kx,ky)=^ { cJ^.^CkCT) vs kx for the ky momenta indi- 
cated, obtained using the ground-state of the two holes 2x8 



cluster and Ck^ 



^E, 



= -»kj , 



where N is the number 



of sites in the cluster, j is a site index and a 
couplings are J=0.2, t'=-0.35 and t"=0.25. 
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The charge structure factor was also calculated in these 
investigations. In Figs.6a-b exact results on a 2x8 clus- 
ter and DMRG results on 2x16 and 2x32 clusters, all at 
hole density x=0.125, are shown. They are very similar, 
illustrating the lack of strong size effects systematically 
observed in our studies of two- leg ladders. The results 
obtained in the more traditionally studied case of J=:0.4, 
t'=t"=0.0 are also shown in Figs.6c-d. They are qualita- 
tively different from those in Figs.6a-b, difference likely 
related with the formation of the quasi-lD regions men- 
tioned before. 
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FIG. 6. Charge structure factor N(kx,ky) vs kx, parametric 
with ky as indicated, (a) corresponds to J=0.2, t'=-0.35, and 
t"=0.25 using the two holes ground-state of the 2x8 cluster 
solved exactly, (b) Same as (a) but using 2x 16 and 2x32 clus- 
ters and 4 and 8 holes, respectively, studied with the DMRG 
technique, (c) Same as (a) but for J=0.4, and t'=t"=0.0. (d) 
Same as (c) but on a 2x16 cluster with 4 holes studied with 
the DMRG technique. Comparing (a) with (b), and (c) with 
(d), it appears that finite size effects are small. 



III. THREE-LEG LADDERS: 

A. One hole: 

Results for three-leg ladders have also been gathered 
in this work. In Fig. 7 the hole is projected at the loca- 
tion shown, from the lowest-energy state at the momenta 
indicated in the caption. In Fig. 7a results are presented 
at J=0.1 and nonzero t'-t". The formation of the across- 
the-hole structure is very clear, in both directions. Along 
the vertical one, a strong AF bond is formed between 
the two spins of the rung where the hole is located. Its 
strength is close to that of a perfect singlet and, thus, 
the coupling of those spins with the rest is small. This 
structure is present at all the couplings investigated as 
shown in Figs.7a-c, even intermediate and large J, and it 
is also clearly present in the exact solution of a three-site 
t-J model at many couplings. In the small J/t regime of 
the standard t-J model, and in the extended t-J model 
as well, this strong AF bond denotes a precursor of the 
TT-shift across the stripe that will form (along the PBC 
direction) as the hole-density increases. Regarding the 
horizontal AF bond (along the legs), its strength depends 
on the particular value of J and t'-t" . In the intermediate 



coupling region often studied, exemplified by J=0.5 and 
no extra hoppings as in Fig. 7c, there is no across-the- 
hole AF correlation along the legs. On the other hand, 
for the couplings of Fig. 7a, the clear ID-like AF segment 
along the central leg induced in the vicinity of the hole 
suggests a larger mobility of the carrier and generation 
of ID-like segments, as discussed in the Introduction and 
for two-leg ladders. 
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FIG. 7. Spin-spin correlations in the state with the low- 
est energy of one hole with momentum k, at the couplings 
described below. The data were obtained exactly on the 
three- leg ladder using the Lanczos method on 3x6 clusters. 
PBC are used in the long direction and open boundary con- 
ditions (OBC) in the short one. Shown are results when the 
hole is projected at the site indicated by an open circle from 
the state under consideration. The dark lines between sites 
represent the absolute value of the (antiferromagnetic) spin 
correlation between the spins located at those sites (scale 
shown). The dashed line in (c) indicates a weak ferromag- 
netic correlation. Results at distance of two lattice spacings 
are shown only near the hole, (a) corresponds to k=(7r/3,0), 
J=0.1, t'=-0.35, and t"=0.25. (b) corresponds to k=(7r/3,7r), 
J=0.1, and t'=t"=0.0. (c) corresponds to k=(7r/3,7r), J=0.5, 
and t'=t"=0.0. 

The tendency toward local spin-charge separation can 
be studied further if the mean- value of the spin in the z- 
direction is calculated when the hole is projected from its 
ground-state to a given site, as done before in Figs.7a-c. 
Since removing one spin from a cluster with an even num- 
ber of sites creates states with spin- 1/2, together with the 
hole in Fig. 8 there must be a spin-1/2 spread over the 
cluster. The results are shown in Fig. 8a for the case of 
small J and nonzero t'-t". In this situation it is clear that 
the z-component of the spin 1/2 is distributed approxi- 
mately uniformly along the central leg where the hole is 
projected (the spin in the outer legs is negligible). In 



agreement with the introductory discussion, the mobile 
hole creates a ID environment (along the central leg in 
this case) to help in its propagation, and in this context 
the spin and charge are separated. On the other hand, 
the results corresponding to an intermediate J of value 
0.5 and no extra hoppings corresponds to a staggered 
spin background surrounding the hole (Fig. 8b), similar 
to results for a vacancy in a Neel background, with the 
exception of the two spins in the same rung as the hole, 
which try to form a strong singlet and thus its mean 
z-component spin tends to vanish. There is a clear qual- 
itative difference between Figs. 8a and b, caused by the 
high mobility of the hole in the presence of a small J 
and nonzero hoppings beyond nearest-neighbors. The 
QP weights Z (not shown) associated to the structure of 
Fig. 8a tend to vanish at k=(0, tt) and others, as found 
before for two- leg clusters and small two-dimensional sys- 
tems. 

(a) _ 




FIG. 8. Mean value of the z-component of the spin ({Sf )) 
denoted by the area of the circles, obtained from the one-hole 
ground-state of the 3x6 cluster with the hole projected at 
the open circle position. Black (gray) circles indicate positive 
(negative) values, (a) corresponds to J=0.1, t'=-0.35, and 
t"=0.25, while (b) corresponds to J=0.5, t'=t"=0.0. 



B. Many holes: 

When two holes are considered on the three-leg ladder 
the situation is similar as in the case already described 
of two-legs, namely each individual hole carries a spin 
arrangement in its vicinity similar to that found in the 
one hole case. Results are shown in Fig. 9a at small J 
and Fig. 9b at intermediate J, in both cases without t'- 
t". In the first case, AF bonds across-the-hole are found 
in both directions, while in the second the bond along 
the legs turns ferromagnetic as for vacancies in a Neel 
background and as in Fig. 7c for one mobile hole |2^. 
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FIG. 9. Same convention as in Fig. 7 but now for the two 
holes ground-state working at: (a) J=0.2, and t'=t"=0.0, and 
(b) J=0.5, and t'=t"=0.0. In this case the most likely con- 
figuration is a tight bound state (not shown), (c) Same as 
(a) but for the three holes ground-state at J=0.25, t'=-0.35, 
t"=0.25. In all cases the momentum is (0,0). 

If three holes are considered at small J (Fig. 9c) the 
most likely hole configuration corresponds to a n^=0.5 
stripe along the central leg, with clear and robust AF 
bonds across it joining the two outer legs (in spite of hav- 
ing open boundary conditions along the rungs). Within 
the central leg or stripe, AF bonds are formed across- 
the-hole, as in a ID system with the same Hamiltonian. 
Fig. 9c qualitatively indicates the dynamical separation 
of a chain subsystem carrying the charge from the rest 
of the spins which are correlated as if that chain would 
not exist. As emphasized in other parts of this paper and 
in Refs. |l^,|l^, this appears to be the way in which the 
system achieves a partial separation of spin and charge 
in two dimensions. 



IV. FOUR- AND SIX-LEG LADDERS: 

A. One Hole: 

Results similar to those found in the case of two- and 
three-leg ladders are also observed with four legs. Re- 
sults for the case of one hole with momentum k=(7r,0) 
are shown in Fig. 10 on a 4x6 cluster solved exactly using 
the Lanczos method. In Fig. 10a results for small J and 



nonzero t'-t" are presented: here the AF bonds across- 
the-hole can be observed and ID-like segments are cre- 
ated near the hole in both directions ||2J]. The effect is 
amplified if at fixed J the values of the hoppings t'-t" are 
increased. Fig. 10b contains the results using abnormally 
large extra hoppings: now the AF bonds across-the-hole 
are quite robust, with a strength that grows as the extra 
hoppings grow in magnitude. There is a smooth con- 
nection between the physically acceptable values of t'-t" 
and those used in Fig. 10b to amplify the effect. Such a 
connection suggests a common origin to the structure. 



(a) 



0.34 

















r^^T^ 









1 1 


. 



PBC 



PBC 



(b) 

















^^?^N 








7^ 





FIG. 10. Spin-spin correlations in the state with the lowest 
energy of one-hole with momentum k, at the couplings de- 
scribed below for a four-leg ladder. The data were obtained 
exactly using the Lanczos method on a 4x6 cluster with PBC 
in both directions. Shown are results when the hole is pro- 
jected at the site indicated by an open circle from the state 
under consideration. The dark lines between sites represent 
the strength of the antiferromagnetic spin correlation between 
the spins located at those sites (scale shown in (a)). Re- 
sults at distance of two lattice spacings are shown only near 
the hole, (a) corresponds to J=0.2, t'=-0.35, t"=0.25 and 
k=(7r,0). Far from the hole the spin correlations are as in un- 
doped ladders, while near the hole they are much distorted. 
Note the clear presence of the across-the-hole antiferromag- 
netic bonds emphasized in this work, (b) Same as (a) but 
for t'=-1.50, t"=1.07 (keeping the same ratio t'/t" as in (a)). 
These large values of the extra hoppings are used to increase 
the magnitude of the effect. 

The one-hole spectral function A(k, uS) is shown in 
Figs.lla-b for the momenta indicated and as a function 
of |t'|. At (0, tt) the QP weight Z rapidly decreases with 
an increasing t' hopping amplitude, while at a momenta 
closer to the expected ground-state momentum of one 
hole it remains more robust, but it eventually tends to 
vanish at large |t'| (see also Fig. 12). These results are 
very similar to those observed for two-leg ladders, for 



planes |12], and also for three-leg ladders (not shown). 



B. Many Holes: 






tVt = - 0l9/\ 


1 ' 1 ' 


3.5'l0''' -~L__ 


f/t = - o.7_y\^ 


8.0 lO" -.^ 


f/t = - 0.5 


A 


1.9 10"',.^ 


f/t = - 0.35 




5.0 10"'^ 


f/t = - 0.2 


^~--— 


3.0 10'^--^ 


f/t = - 0.1 




7.9 10"'~-„ 

A, 


f/t = 0.0 

1 


1 . J r- 


1.3 10-'--. (a) 

— (^ . .1 — . — Ji\ . 




FIG. 11. One-hole spectral function A(k,(i;) obtained ex- 
actly from the undoped ground state of the 4x6 cluster with 
PBC in both directions, working at J=0.2, with t'=-1.4t" 
fixed, and parametric with t'. The arrows indicate the po- 
sition of the first pole in the spectra and its weight Z (nor- 
malized such that 0<Z<1). (a) corresponds to k=(0,7r) and 
(b) to k=(7r/3,7r/2) (the closest to k=(7r/2,7r/2) in the cluster 
considered here). 
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FIG. 12. Weight Z of the lowest energy state in the A(k, lj) 
one- hole spectral function of exactly solved 4x6 clusters at 
J=0.2, t'=-1.4t", as a function of |t'/t|. Z is normalized such 
that 0<Z<1. Shown are results at the momenta indicated. 



To clarify the behavior of two holes on four-leg ladders, 
consider the use of cylindrical boundary conditions (OBC 
in one direction and PBC in the other [|j). Suppose 
the two-hole ground-state is considered and one hole is 
projected to a site belonging to the central legs, where 
the density of holes is the largest. In this situation it is 
possible to study the probability of finding the other hole, 
with results shown in Fig. 13a: the second hole clearly 
prefers to be along the PBC direction and still within 
the two central legs. Actually the largest probability is 
found at a distance of two lattice spacings along the same 
leg. The results can be interpreted as the formation of 
a loop of charge which wraps around the direction with 
PBC, as discussed recently in Ref. [^ and also in Ref. 
Q. The state appears to correspond to a n/i=0.5 stripe, 
not rigid but fluctuating in the direction perpendicular 
to it. 
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FIG. 13. Results corresponding to the exact two holes 
ground-state of a 4x4 cluster with PBC in one direction and 
OBC in the other, as indicated. The couplings are J=0.2, 
t'=-0.35, and t"=0.25. (a) corresponds to the case where one 
hole is projected from the ground-state into the location de- 
noted by the open circle. The full circles have an area propor- 
tional to the probability of finding the second hole at a given 
site. A loop of mobile holes is dynamically formed around the 
PBC direction, (b) AF spin-spin correlations with the con- 
vention followed in previous figures, for the case where the 
two holes are projected from the ground-state at the location 
indicated by open circles. This hole configuration is the one 
with the largest weight in the ground-state. Results at dis- 
tance of two lattice spacings are shown only near the stripe. 

In Fig. 13b the spin correlations are shown for the 
holes projected into the configuration with the largest 
weight in the two-hole ground-state. In agreement with 
the discussion associated with Fig. 13a, the holes appear 
to be forming a n;i=0.5 stripe along the direction with 
PBC. Moreover, the spin correlations across-the-stripe 
are clearly antiferromagnetic as in the three-leg ladder 
(Fig. 9c) and in experiments 0], and the two spins of the 
four-site stripe considered here are correlated also anti- 
ferromagnetically. This result is representative of a large 
number of similar data gathered in our present analy- 
sis, namely it appears that in general holes tend to form 



loops of charge around the closed direction if a system has 
cylindrical boundary conditions. This is also in excellent 
agreement with the conclusions of earlier work by White 
and Scalapino ^, although their description of stripes 
is sometimes based on the condensation of d-wave pairs 
|[7| while ours is based on a kinetic energy optimization 
(namely a one-hole problem). 

The appearance of stripes based on the 4x4 results of 
Fig. 13 needs to be confirmed increasing the cluster size. 
This analysis was done in part in Ref. |13| as discussed in 
the Introduction, but here those results are expanded and 
more details are provided. Note that our present analysis 
using the DMRG method is restricted to the t-J-t' model, 
namely the t" hopping will be considered to be zero. The 
reason is that in the implementation of the DMRG tech- 
nique sites are aligned along a one-dimensional pattern 
even for ladder geometries, and a t" hopping would link 
sites along this equivalent chain which are located several 
lattice spacings from each other, reducing the accuracy of 
the method. Nevertheless even without t" the stripes re- 
ported in Fig. 13 and Ref. Q can be clearly observed on 
larger systems. For example, consider in Fig. 14 two holes 
on a 4x12 cluster with cylindrical boundary conditions 
(PBC along the rungs), J=0.5, and t'=-0.3 (t"=0.0). One 
of the holes is projected to an arbitrary site of one of the 
central rungs, where the hole density is the largest. The 
distribution of the second hole around the projected one 
is quite similar to the result observed on the 4x4 clus- 
ter, namely the largest chances are at distance of two 
lattice spacings along the rung. This state resembles a 
bound-state in the sense that the two holes are close to 
each other, but its shape is better described as a stripe 
configuration or a loop of charge that wraps around the 
short direction. 



a 4x8 cluster at small J and without t' and t". In Fig. 15 
the density of holes is shown for the case where one hole 
is projected into one of the rungs with the largest den- 
sity. It is clear from the figure that in the vicinity of 
the projected hole there is a structure similar to that 
observed for the two-hole ground state, namely a local 
maximum in the hole density is observed at two lattice 
spacings from the projected hole. In addition, clearly a 
large accumulation of charge appears on the other side 
of the cluster from where the projected hole is located. 
This other sector populated by holes involves two rungs 
in width, and it is framed in Fig. 15. Note that the den- 
sity in that second stripe is very uniform showing that 
there are no charge correlations between the two stripes 
(the two charged loops are independent of each other). 
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FIG. 15. Results corresponding to the four holes 
ground-state treated with the DMRG technique on a 4x8 
cluster with PBC in the short direction and OBC in the other, 
as indicated. The couplings are J=0.2, and t'=t"=0.0. One 
of the holes is projected at the open circle site. The full cir- 
cles denote the probability of finding another hole at a given 
site. Two loops of mobile holes appear to be formed wrapping 
around the PBC direction (one of the two is framed) . 
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FIG. 14. Results corresponding to the two holes 
ground-state studied with the DMRG technique on a 4x12 
cluster with PBC in the short direction and OBC in the other, 
as indicated. The couplings are J=0.5, t'=-0.3, and t"=0.0. 
One of the holes is projected at the open circle site. The full 
circles denote the probability of finding the second hole at a 
given site. A loop of mobile holes is formed wrapping around 
the rung direction. 

If more holes are added to the four-leg ladder, it ap- 
pears that stripes similar to those observed in Fig. 14 are 
formed. For example, consider the case of four holes on 



In Fig. 16a the charge structure factor is shown for the 
case of an 4x8 cluster with four holes, studied with the 
DMRG method at small J. The result resembles those 
obtained for two-leg ladders in Fig. 6. In Fig. 16b the spin 
structure factor is shown for the same cluster and density. 
Spin incommensurability clearly appears in this system, 
as remarked before in Ref. p3| working at other cou- 
plings, and as presented for the two-leg ladder (Fig. 5a) 
as well. It is clear that the two-, three- and four-leg lad- 
der systems share very similar physics in the charge and 
spin sectors. 
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FIG. 16. Results obtained with the DMRG technique on 
a 4x8 cluster with four holes (x=0.125) at J=0.2, and 
t'=t"=0.0, using cylindrical boundary conditions (PBC along 
rungs), (a) Charge structure factor N(kx,ky) vs kx, paramet- 
ric with ky as indicated, (b) Spin structure factor S(kx,ky) 
vs kx, parametric with ky as indicated. Incommensurability 
is clear in the ky=7r branch. 



C. Six-leg Ladders: 

Accurate results for six-leg ladders are difficult to ob- 
tain in the regime of parameters studied in this paper, 
namely small J and nonzero t'-t". However, some non- 
trivial results can still be gathered. For instance, consider 
in Fig. 17 the case of three- holes on a 6x4 cluster with 
PBC along the direction with six sites, and OBC along 
the short one. In this case a good approximation to the 
ground-state can be obtained using the ORBA method 
with about one million 4-site rung-basis states. The pro- 
cedure to obtain the results shown in Fig. 17 are the fol- 
lowing: first it was noticed that the central long legs are 
the ones the most populated by holes. Then, one hole was 
projected at an arbitrary site belonging to those central 
legs. In the next step, the hole density is obtained with 
that hole projected, and at the position with the largest 
density a second hole is now projected. The distance be- 
tween the first and second projected holes is two lattice 
spacings. With those holes projected as shown in Fig. 17, 
the hole density is recalculated. It is clear that at two 
lattice spacings from the projected holes the density has 
a maximum. This result is once again compatible with a 
Uh=0-5 stripe formed this time by three-holes in a closed 
loop around the direction with PBC. The stripe is not 
rigid but fluctuating perpendicular to its main direction. 
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FIG. 17. Results corresponding to the three holes 
ground-state treated with the ORBA technique on a 6x4 clus- 
ter with PBC in the long direction and OBC in the other, as 
indicated. The couplings are J=0.2, t'=-0.35, and t"=0.25. 
Two of the holes are projected at the open circle sites. The 
full circles denote the probability of finding the third hole at 
a given site. A loop of mobile holes appears to form wrapping 
around the PBC direction. 

The density of holes with a given momentum k for the 
6x4 cluster with three holes is shown in Fig. 18. Clearly 
the momenta the most important are those around (tt, tt), 
result similar to those reported for the two-legs cluster in 
Fig. 5b and for a 4x6 two-holes cluster with PBC in both 
directions in Ref. ilSl . 




FIG. 18. n(kx,ky)=X;<.( cLck^) vs kx for the ky 

momenta indicated, obtained using the ORBA approxima- 
tion to the ground-state of three holes on a 6x4 cluster 
with PBC in the x-direction and OBC in the y-direction. 
'*''■' Cio-, where N is the number of sites in 
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the cluster, j is a site index and a = +, — . 
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V. ANALYTICAL CALCULATIONS: 
SPIN-POLARON APPROACH 

In previous sections, it has been shown among other 
items that the across-the-hole correlations identified be- 
fore numerically on small square clusters for the one-hole 
problem ||l2[ exist also on ladders. For completeness, in 
this section an analytical approach to the one-hole prob- 
lem in two dimensions is described. In agreement with 
the above mentioned results, once again the existence of 
across-the-hole correlations is confirmed this time using 
a non-numerical method, highlighting the robustness of 
this feature. 

The analytic approach used here is based on the picture 
of the spin-bag QP or magnetic polaron p^|2^ . The spin 
deviations from a perfect antiferromagnet in the vicinity 
of the hole are described by a set of fiuctuation opera- 
tors An | |2^ , |27| , p8[ | , which create strings of spin defects 
attached to the hole. The one-hole Green's function is 
evaluated using a cumulant version p9[ | of Mori-Zwanzig 
projection technique. 

The undoped ground state is modeled by 1-0) = 
exp(^^ Q;yS'y)|0Neei) where the operators 5^ create clus- 
ters of spin fluctuations in the classical Neel state. Fol- 
lowing Rcf. iQ it is possible to obtain a non-linear set 
of equations for the coefficients a^, which is solved self- 
consistently. The fluctuation operators An for describ- 
ing the spin deviations caused by hole motion include 
the usual string operators up to a certain length ^max 
and additional operators for local fluctuation configura- 
tions. The Green's function is calculated by projection 
technique in the subspace spanned by the operators A„, 
which amounts to the diagonalization of the dynamical 
matrix (tplAlj^HAnlip). Processes outside the subspace 
formed by An\ip) (i.e., self-energies) are neglected, there- 
fore a discrete set of poles for each spectrum is obtained. 
The eigenvectors of the dynamical matrix represent the 
wavefunctions corresponding to the poles and they al- 
low for the calculation of expectation values such as spin 
correlation functions. In the following the pole next to 
the Fernu surface representing the QP is the only one 
considered, the remaining part of the spectrum forms an 
incoherent background well separated from the QP peak. 
The approach sketched here has been successf ully ap plied 
to the one-hole problem in several contexts ||2^,^; de- 
tails are described in Ref. pq] . In the present calculations 
strings with /,nax — 5 have been used (total number of 
fluctuation operators An around 1200). Note that the 
accuracy of the approximation decreases with increasing 
t since the employed expansion parameter is basically 
the amplitude of the deviations from the undoped state. 
Therefore it is better to restrict ourselves to the param- 
eter range with J/t>0.05 where the comparison with nu- 
merical results reveals the good accuracy of the analytical 
approach. 

First, let us discuss the weight of the QP pole which 
is known to decrease with decreasing J/t since the size 



of the spin polaron increases. The full parameter depen- 
dence of the QP weight Z at momentum (tt, 0) is shown 
in Fig. 19a. The introduction of t'<0 rapidly suppresses 
the value of Z^o even for intermediate or large J, such 
that the entire region of moderate t' is characterized by 
small Ztt^o, in good agreement with the computational re- 
sults shown in this paper and in Ref. |l3] . On the other 
hand, the t' dependence of Z at (tt/2, 7r/2) is weaker (not 
shown), also as found in Ref. [Of. Nevertheless, its asso- 
ciated QP weight Z is also strongly suppressed with the 
reduction of J. 
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FIG. 19. (a) QP weight Z and (b) magnitude of 
across-the-hole correlations Sah at momentum (7r,0) in the 
plane J/t vs. t'/t, with t'/t" fixed to -1.4. 

In order to study the spin configuration near the mo- 
bile hole, static spin correlations in real space have been 
calculated relative to the hole position. It is known 
that the limit of small t (static vacancy) leads to an 
increase of the antiferromagnetic correlations on bonds 
near the hole, whereas hole hopping tends to "scram- 
ble" the spins in the hole environment which in general 
weakens the AF tendencies. This general behavior is 
well reproduced by our calculations. However, the in- 
troduction of t' leads additionally to antiferromagnetic 
correlations Sah across-the-hole, see Figs. 19b and 20b, 
quite consistent with those found numerically here and 
in previous studies. In Fig. 20b, the strength of the an- 
tiferromagnetic bonds at J=0.05, t'=-0.35, t"=0.25, and 
momentum (tt, 0) are shown as illustration. Across the 
mobile hole strong antiferromagnetic correlations develop 
which are supported by further strong AF bonds form- 
ing a chain segment. The other nearest-neighbor bonds 
are weaker by a factor of two than the chain bonds. The 
values of Sah , obtained with the analytic approximation 
used in this section are shown in Fig. 19b. Although they 
are not as strong as found numerically, the qualitative 
trends agree quite well with computational studies. 

The stronger tendency for local spin-charge separation, 
at small J/t or in the presence of t' and t", can be un- 
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derstood within the string picture. The string of defect 
spins attached to the hole basically connects the spin 
and charge parts of the excitation, confining them at 
long distances. Therefore longer average strings corre- 
spond to less tightly bound spin and charge, i. e., larger 
spin polarons, immediately implying a smaller quasipar- 
ticle weight. The energy cost per unit length of string 
is proportional to J, hence smaller J/t allows for longer 
strings [Q . Interestingly, the presence of t' and t" leads 
to a related effect, namely the hole can easily 'get rid' of 
its attached string by intra-sublattice hopping, the string 
being then absorbed in the background spin fluctuations. 
This mechanism makes long strings effectively less costly, 
leading to the described large polarons with small QP 
weight Z. 
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FIG. 20. (a) Ratio of the effective QP masses at mo- 
mentum (±7r/2, ±7r/2) which is a direct measure of the 
(an)isotropy of the dispersion. In the region near t'=0 and 
J/t=0.4 the value is around 0.14 (in agreement with previous 
numerical studies [16]). (b) Illustration of the static spin cor- 
relations near the hole at J=0.05, t'=-0.35, t"=0.25 and mo- 
mentum (tt, 0). The width of the lines is proportional to the 
spin-spin correlation between sites i and j, namely —(Si ■ Sj), 
and the value of the across-the-hole correlation is -0.1. 

For completeness, in Fig. 21 the calculated QP disper- 
sions along high-symmetry lines in the Brillouin zone are 
shown. The points correspond to (a) J=0.40, t'=0.0, (b) 
J=0.40, t'=-0.35 and (c) J=0.05, t'=-0.35 with t'/t''^- 
1.4. The change in the curvature near (±7r/2, ±7r/2) is 
clearly visible when adding t' and decreasing J. Whereas 
the dispersion is strongly anisotropic for large J and small 
t', there is an entire region of near isotropy for moder- 
ate values of t' and a large range of J/t, in agreement 
with PES experiments 11^1. This is also illustrated in 
Fig. 20a where the ratio of the effective QP masses at 
Ko=(7r/2,7r/2) for the entire parameter space are pre- 
sented. The masses are defined as usual as eigenvalues of 
the tensor m given by ek = (2m);;^^(k - Kq) a{\<i - 'Kq) p 
where et is the QP energy near Kq. 




(0.0) (1,1) (7t,7t) 



(0,7t) (0,0)(0,7t) (|,|) (71,0) 

Momentum 



FIG. 21. QP dispersions for points (a), (b) and (c) as men- 
tioned in the text. The energy zero level has been set at the 
center of mass of the band. The evolution towards a nearly 
isotropic maximum around (±7r/2, ±7r/2) is clearly visible. 

Summarizing, the introduction of t' leads to a suppres- 
sion of the QP weight, especially at momentum (tt, 0), a 
delocalization of the spin carried by the spin polaron, and 
the formation of antiferromagnetic correlations across the 
mobile hole. All the results are in good qualitative agree- 
ment with the numerical calculations reported here. 



VI. CONCLUSIONS 

In this paper, tendencies toward spin-charge separation 
at short distances have been discussed in the context of 
the extended t-J model using ladder geometries. This ef- 
fort generalized previous calculations carried out on small 
square clusters. Analytic approximations have also been 
used, with results in good agreement with the computa- 
tional ones. Overall it is concluded that in regimes where 
the hole kinetic part of the Hamiltonian dominates, holes 
tend to arrange the spin environment in such a way that 
across-the-hole robust antiferromagnetic correlations are 
generated, both on ladders and planes. This arrangement 
helps the hole move easily among the spins, and it resem- 
bles the structure found in one-dimensional spin-charge 
separated systems. For a variety of reasons described 
here, it is believed that at least at short distances similar 
tendencies toward spin-charge separation are at work in 
ladders and two-dimensional systems, at very small J in 
the standard t-J model, or in the extended t-J model. 
At finite hole density, holes share their nontrivial spin 
environment, forming half-doped stripes as recently dis- 
cussed by some of the authors in Ref . ||l^ . Here more ev- 
idence substantiating this previous result has been pro- 
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vided. The stabilization of stripe tendencies discussed 
here is based upon a small J/t picture, and it has no ob- 
vious relation with those emerging in the opposite limit 
of large J/t based upon the frustration of phase sepa- 
rated tendencies. Although more work is certainly still 
needed to confirm the picture described here, it appears 
that a gas of spinous and holons, as envisioned in two- 
dimensional theories for cuprates based upon spin-charge 
separation, may not be operative at low hole densities, 
but instead stripes of holons seem to be forming. 
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